First Surface-resolved Results with the IOTA Imaging 
Interferometer: Detection of Asymmetries in AGB stars 



Ragland, S. 1 ' 2 , Traub, W.A. 3 ' 1 , Berger, J.-P. 4 , Danchi, W.C. 9 , Monnier, J. D. 6 , Willson, L. 
A. 11 , Carleton, N. P. 1 , Lacasse, M. G. 1 , Millan-Gabet, R. 5 , Pedretti, E. 6 , Schloerb, F. P. 7 , 
Cotton, W. D. 8 , Townes, C.H. 10 , Brewer, M. 7 , Haguenauer, P. 12 , Kern, P. 4 , Labeye, P. 13 , 
Malbet, F. 4 , Malin, D. 7 , Pearlman, M. 1 , Perraut, K. 7 , Souccar, K. 7 , Wallace, G. 7 

ABSTRACT 

We have measured non-zero closure phases for about 29% of our sample of 
56 nearby Asymptotic Giant Branch (AGB) stars, using the 3-telescope Infrared 
Optical Telescope Array (IOTA) interferometer at near-infrared wavelengths (H- 
band) and with angular resolutions in the range 5-10 milliarcseconds. These non- 
zero closure phases can only be generated by asymmetric brightness distributions 
of the target stars or their surroundings. We discuss how these results were 
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obtained, and how they might be interpreted in terms of structures on or near 
the target stars. We also report measured angular sizes and hypothesize that 
most Mira stars would show detectable asymmetry if observed with adequate 
angular resolution. 

Subject headings: stars: AGB, stars: asymmetric stars, stars: circumstellar 
shell, stars: surface features, stars: non-radial pulsation, technique: long baseline 
interferometry, technique: closure phase 

1. Introduction 

The stars in this study are all AGB stars, that is, stars found at or near the tip of the 
Asymptotic Giant Branch in the HR diagram. They are low to intermediate mass stars, 
having already spent most of their lives as normal stars, and currently heading towards their 
deaths probably in the form of planetary nebulae, leaving the central star as a white dwarf. 
Most AGB stars are variable in brightness; those with relatively regular and large amplitude 
visual variations (> 2.5 mag) with periods in the range 100 - 1000 days are classified as Mira 
variables. The Miras and some of the other, semi-regular (SR) or irregular (Irr) variables 
have observed mass loss rates ranging from 10~ 7 to > 10 -5 M Q per year (Knapp & Morris 
1985). Diameter changes, opacity changes, and possibly other processes such as convection 
contribute to the brightness variation in these stars. 

The Mira stage of evolution has been identified as marking the onset of the "superwind" 
phase, i.e. that evolutionary stage where mass loss rates rapidly increase and result in the 
termination of AGB evolution (Bowen & Willson 1991, Willson 2000). These stars thus 
serve as markers for the tip of the AGB in various populations, something already known 
for the shorter period cases from the few Miras that appear in globular clusters such as 47 
Tuc (Frogel, Persson & Cohen 1981). 

Miras with close companion white dwarfs usually are classified as symbiotic systems 
(Allen 1984; Whitelock 1987; Luthardt 1992; Belczyhski et al. 2000). A few Miras are known 
to have companions but are not (or are only very mildly) symbiotic systems; this includes o 
Cet = Mira, with a probable WD companion in a multi-century orbit (Reimers & Cassatella 
1985; Wood & Karovska 2004). Statistics for the binarity of Miras are otherwise quite 
uncertain, in part because the expected orbital velocity amplitudes for a close companion, 
30 km/s at 1 AU and 10 km/s around 5 AU, are very similar to the shock amplitudes of 
20-30 km/s produced by the Mira pulsation itself (Hinkle, Scharlach & Hall 1984). 

In this paper we use the word "asymmetry" to mean that part of the 2-dimensional 



- 3- 



brightness distribution which cannot be made symmetric with respect to a reflection through 
a point. Thus, for example, an elliptical uniform disc or an equal-brightness binary system 
are both symmetric, but a binary with unequal brightness or a star with an off-centered 
bright/dark spot is asymmetric. 

Departure from circular symmetry has been known in AGB stars from various high 
angular resolution observations (Karovska et al. 1991; Wilson et al. 1992; Haniff et al. 1992; 
Richichi et al. 1995; Ragland 1996; Weigelt et al. 1996; Karovska et al. 1997; Tuthill et al. 
1997; Lattanzi et al. 1997; Wittkowski et al. 1998; Tuthill et al. 1999, 2000; Hofmann et al. 
2000; Thompson et al. 2002; Monnier et al. 2004a; Weiner et al. 2006). The observed depar- 
tures from circular symmetry have been interpreted either in terms of elliptical distortions 
or an otherwise symmetric photosphere containing localized compact features. However, 
no consensus exists as to the mechanism that would cause such departures from apparent 
circular symmetry. 

Dust shells surrounding AGB stars have observed asymmetries as well such as in Mira 
(Lopez et al. 1997), the carbon stars IRC+ 10216 (CW Leo; Tuthill et al. (2000)), CIT 6 
(Monnier, Tuthill & Danchi 2000) and IK Tau (Weiner et al. 2006) among others. The 
connection between apparent surface features and the morphology of the dust shells has not 
been established. 

About 50% of all planetary nebulae (PN) display bipolar symmetry (Zuckerman & Aller 
1986), but only a small fraction of circumstellar envelopes show bipolarity. A surprisingly 
large number of proto-PN show roughly circular arcs surrounding a bipolar core, suggesting 
that in most cases the AGB mass loss is spherically symmetric and the asymmetry seen in 
the PN occurs well after the Mira stage (Su 2004; Willson & Kim 2004). Recent studies of 
jets around a few AGB stars (Kellogg, Pedelty & Lyon 2001; Imai et al. 2002; Sahai et al. 
2003; Sokoloski & Kenyon 2003; Brocksopp et al. 2004) from radio, x-ray or Hubble Space 
Telescope (HST) observations suggest that those stars showing substantial asymmetry may 
all have a low mass stellar companion accreting mass from the AGB primary. Recent SiO 
maser observations of AGB stars show departures from spherical symmetry (Diamond et al. 
(1994), Greenhill et al. (1995), Diamond & Kemball (2003), Cotton et al. (2004), Soria-Ruiz 
et al. (2004)). The observed circumstellar SiO masers tend to occur in clumpy, partial rings 
centered on the central star (Diamond et al. (1994)). Cotton et al. (2002) observed 9 stars 
in SiO, at least two of them known binaries, and used the modeling of Humphreys et al. 
(2002) in discussing the results. SiO maser emission comes from ~ 2AU or ~ 2R* where 
the outflow velocity gradient along the line of sight is small. The special conditions required 
for maser emission potentially gives rise to bias in the statistics of asymmetry in the sample 
population of stars. 
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In this paper we report the initial results from one phase of a larger program, the 
Mira Imaging Project, to investigate asymmetries in AGB stars using three interferometer 
facilities, each capable of making closure phase measurements. These facilities are the IOTA, 
the Infrared Spatial Interferometer (ISI), and the Very Long Baseline Array (VLBA). The 
present paper focuses on IOTA results. 

Subsequent to the work reported here, and as a part of the ongoing Mira Imaging 
Project, selected Mira targets with positive closure phase signal from our survey have been 
re-visited at different pulsational phases, baselines, position angles, and wavelengths in order 
to characterize the observed asymmetry. The results of this ongoing study will be presented 
elsewhere. In this article, we present the initial survey results for all our targets. 

2. Observations 

The observations reported here were carried out during the commissioning phase of the 
IOTA 3-telescope array (Traub et al. 2004) and integrated-optics beam-combiner, IONIC 
(Berger et al. 2004), operating in the H-band (1.65 /zm) atmospheric window. Observations 
of binary stars taken with the same instrumental configurations were reported by Monnier 
et al. (2004b) and Kraus & Schloerb (2004). We report here the results of the first phase of 
our program in which we have studied 56 evolved giants (Tables 1 & 2) including 35 Mira 
stars, 18 SR variables and 3 Irr variables 1 looking for asymmetry in their brightness profiles. 

We report observations taken during six observing runs during May 2002 to May 2003. 
Observations were taken either with a standard H band filter (A G = 1.65 /xm, AA = 0.3 
(Ltm) or with a narrow-band filter (A Q = 1.64 /im, AA = 0.1 fim). Typically, five minutes of 
program star observations were followed by nearby calibrator observations under identical 
instrumental configurations. For the observations taken during March 2003 and May 2003 
observing runs, we used ND filters for the bright targets since we had excellent optical 
throughput with newly coated primary mirrors and well optimized beam-train. On each 
star, we record 4 sets of data files each containing about 500 scans. A scan consists of 
changing the optical path difference between two beams by roughly 75 /im in saw tooth 
form. We then take about 400 scans of shutter data for calibration. The shutter data 
sequence consists of allowing only one beam at a time (telescope A, B and then C) and at 
the end blocking all three beams. Each scan takes about 100 ms. 



1 82 Lyr is classified in the Combined General Catalogue of Variable Stars (CGCVS; Samus et al. (2004)) 
with an uncertainty as a SR variable, and no period estimation is available in the literature. We consider 
this target as being an Irr variable for the purpose of this paper. 
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All targets were observed with a three-baseline interferometer configuration, forming 
a closed triangle. Earth rotation enables closure-phase measurements at slightly different 
projected baselines (and hence different closed triangles) when observations are made at 
different hour angles. We have adopted baseline bootstrapping (Mozurkewich & Armstrong 
f 992) at the IOTA whereby fringes are tracked on two short baselines, while science data are 
recorded on all three baselines simultaneously, enabling low visibility measurements on the 
third (long) baseline. Details of the detector camera and the fringe tracker algorithm used 
for this work are reported by Pedretti et al. (2004) and Pedretti et al. (2005) respectively. 

IOTA's maximum baseline of B=38 m yields an angular resolution of A/2B ~ 4 mas at 
1.65 /iin. The present faint limit with the IONIC beam-combiner is H ~ 7 for the broad band 
filter and H ~ 5 for the three narrow band filters. For the difficult case of observing well 
resolved Mira stars with the H filter at or below 5% visibility level, the limiting magnitude 
is H ~ 4. The limiting magnitude of the star tracker at IOTA, for these observations, was 
V ~ 12 for late-type stars. The angular resolution of our short south-east arm at IOTA is 
lower than that of the North-east arm, meaning that we could miss some asymmetry if it 
were predominantly parallel to the projected south-east baseline of the interferometer. 

Pulsation periods for all Mira and SR variables (53 out of 56 program stars) are from 
Combined General Catalog of Variable Stars (CGCVS). Among these 53 program stars, 37 
also have period estimations from American Association of Variable Star Observers (AAVSO) 
derived using a data window centered on JD 2452000 (Templeton, Mattei & Willson 2004) 
enabling us to validate the CGCVS periods. In addition, AAVSO has tentative or very 
tentative periods for three more SR stars. The CGCVS periods are consistent with available 
AAVSO periods for all but four of our program stars. Interestingly, the AAVSO periods for 
all four discrepant stars, namely X Cnc, BG Ser, UU Aur and W Ori, are roughly twice that 
of CGCVS periods although the AAVSO periods for two of them, namely UU Aur and W 
Ori, are either tentative or very tentative values. For completeness, one of the Irr variable, 
namely, TX PSc has a very tentative AAVSO period of 255.5 days. 



3. Data Reduction 

The recorded interferograms were reduced with an IDL code package developed by one 
of us (Ragland). Our single- mode integrated-optics beam-combiner chip (Berger et al. 2004) 
has 3 input beams (I a , h, I c ), and 6 output beams (ij, i — 1 — 6). Each input is split 
into 2 parts and coupled to the outputs as follows: (a, b) <-> (1,2); (a, c) <-> (3,4); (b,c) <-> 
(5,6). The complimentary outputs (I2, I4 & I&) have the same information as the normal 
outputs (Ji, I 3 & J 5 ) except for a n fringe intensity phase shift. Hence the normal and the 
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complimentary outputs could be combined in order to improve the signal to noise ratio of the 
measurements. The background subtracted outputs are combined two by two (with opposite 
signs) and normalized as follows, in order to remove scintillation noise which is common to 
both normal and complimentary outputs, 



t h/!l ~ h,h m 

lab ~ ~ , (l) 

and correspondingly for the other outputs. Here I denotes the mean over the entire scan. 

The power spectra of the resultant three outputs I ab , I bc and I ac are computed and the 
fringe power for each of these three outputs is estimated by integrating the power (P) un- 
der the fringe profile after background power subtraction (Baldwin et al. 1996). The fringe 
power is proportional to the visibility-squared (V 2 ). The target V target * s calibrated by mea- 
suring the fringe power for a nearby calibrator of known V 2 aH5 under the same instrumental 
configuration and by taking the ratio, i.e. 



V 2 = V 2 [ (o) 

target calib \ p J v > 

\ * calib J 

The closure phase is the sum of the fringe phases simultaneously observed on three 
baselines forming a closed triangle and is insensitive to phase errors induced by the turbulent 
atmosphere or optics (Jennison 1958). If the phase errors introduced into the three beams 
are S a , 5 b and 8 C , then the observed fringe phase between baselines a and b can be written as 

(Pab = 4>ab + 8 b -5 a (3) 

Here ip ab is the true object fringe phase between baselines a and b. 

The observed closure phase $ c / is equal to the true object closure phase, i/j a b + ?Pbc + ?Pca, 
to within the measurement noise, as follows: 



$d = <\>ab + <Pbe + <j>ea + noise (4) 

= ipab + S b - 5 a + ipbc + 5 C - 5 b + ip ca + 5 a - 5 C + noise (5) 
= ^ab + ipbc + ipca + noise (6) 

We estimate closure phase as the phase of the bispectrum (Weigelt 1977) of simultaneous 
fringes obtained with the three baselines. The instrumental closure phase is estimated using a 
nearby calibrator and subtracted from the raw closure phase of the target to give a calibrated 
target closure phase. 
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Typical one-sigma formal errors in our uncalibrated closure phase and V 2 measurements 
are ~ 0.2° and ~ 2% respectively. The formal errors are estimated from the scatters of the 
500 fringe scans. In the case of V 2 errors, the error due to background power subtraction 
is also incorporated into the formal error. The calibration process adds up additional errors 
and the measurement error is estimated as the square-root of the sum of the squares of the 
formal and calibration errors. The one-sigma measurement error is reported in Tables 1 & 
2. We estimate calibration errors by observing calibrators under same observing condition 
and calibrating one calibrator with the other after accounting for the finite sizes of both 
calibrators. The estimated one-sigma calibration error is 5% for the V 2 measurements and 
0.5° for closure phase measurements. There could be unaccounted systematic errors in our 
measurements. For the purpose of this article, we adopt a systematic error of 2° for our 
closure phase measurements. The total error is estimated as the square-root of the sum of 
the squares of the formal, calibration and systematic errors. If the measured closure phase is 
less than twice the total error in the measurement then we call it essentially a non-detection 
of asymmetry. However, if it is larger, we call it positive detection. Further discussion on 
our visibility and closure phase measurements could be found in Ragland et al. (2004). 



4. Results 

Targets with centro-symmetry should give a closure phase of either zero or ±180° de- 
pending on how many baselines are beyond the 1st, 2nd, etc nulls. The majority of our 
targets show zero closure phase. However, 14 of the 56 have non-zero closure phase. Of 
these, 12 are Mira stars, and 4 are SR/Irr variable carbon stars. Among the 12 Mira stars, 
all but x Cyg are oxygen-rich Mira stars; x Cyg is classified as an S star. The frequency of 
asymmetry from our studies is 34% in Mira stars, 17% in SR variables, 33% in Irr variables 
and thus 29% in our entire sample of AGB stars. In terms of chemistry, the frequency of 
asymmetry is 33% in carbon stars and 27% in oxygen-rich stars. 

Table 1 gives the measured closure phases along with observational and target informa- 
tion for these 14 targets that show measurable asymmetry from our observations. We have 
included in this table one observation each for R Cnc, R LMi & V Hya in which no asym- 
metry is detected since we use these observations for size estimation. All targets (except 
R Aur, x Cyg & UU Aur) have follow up measurements taken with all three narrow band 
filters to characterize the observed asymmetry. These results will be presented elsewhere. 
The V 2 data for these targets are fitted with uniform disk (UD) models (figure 1) and the 
measured angular sizes are also given in Table 1. Table 2 lists targets for which we have 
not detected asymmetry from our survey. The calibrators used for our observations are also 
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listed in these tables. The angular sizes for most of the calibrators are taken from Wesselink, 
Paranya & de Vorkin (1972). Typically, we use calibrators with angular sizes less than 3 
mas for our measurements. However, during the early part of the survey, we had to use 
larger calibrators because of low throughput of the instrument. We adopted the measured 
sizes from interferometric technique for a Vul and 7 Sge (Hutter et al. 1989), and from lunar 
occultation measurements for UU Aur (Bohme 1978). We estimated the angular size of 7 
Peg from V and K magnitudes (van Belle 1999). 



Table 1. Derived closure phase and UD diameters in the H band for targets with detected asymmetry. 
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Table 1 — Continued 
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3234.3 
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^his observation doesn't show asymmetry. 

2 Observed with the narrow band filter at 1.64//m. 

3 At least one of the baseline resolved this target to the level of below 5% in V 2 . 



4 UD model failed to fit the data. The derived value probably gives the diameter of the dust shell rather than diameter of 
central star. 

5 The uncertainties indicates random errors only (see text). 
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Six of our targets with detected asymmetry, namely, U Ori, R Cnc, R LMi, S CrB, R 
Aql, R Aqr have earlier H band size measurements (Millan-Gabet et al. 2005) taken within 
± 0.2 pulsation phase with respect to the pulsation phase of our measurements for these 
stars. We have plotted our size measurements against the measurements by these authors 
for these six targets in Figure 2. The scatter in this figure is comparable to the scatter among 
measurements of a given star at multiple epochs in a single program, and thus probably either 
signifies actual size variation at the source or that the interpretation of the observations (UD) 
is too simple. Several targets have earlier K band angular size measurements (van Belle et 
al. (1996); van Belle, Thompson & Creeck-Eakman (2002); Millan-Gabet et al. (2005); 
Mennesson et al. (2002); Dyck, van Belle & Benson (1996)). 
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Fig. 1. — Visibility-data fitted with an UD model for the targets with positive asymmetry detection from 
our observations. UD models fits the data very well except for V Hya. In the case of V Hya, the derived size 
is possibly the size of the dust shell rather than that of the central star. 
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Shows the comparison of sizes reported in this article with those reported in the literature. 



Table 2. Closure phase measurements for targets with no detectable asymmetries 



CGCVS AAVSO 



Target 


Date 


Calib. 


Period 


Period 


Phase 


Sd Tv 


D max 




Yep 




UT 




( davs") 








( niasl 


(dee) 




(a) Mira Stars: 


















U Per 

W J. CI 


30Oct02 


51 And 


320.3 


318.7 


0.5 


M6e 


3607.2 


5.84 2 


1.4±1.1 


R Tri 


29Tan03 


24 Per 


266.9 


264.8 


0.9 


M4IIIe 


3548.7 


4.46 3 


2.8±0.5 


RT Eri 


31 Tan03 


HR 1543 


370.8 


376.5 


0.9 


M:e 


2527.9 


6.3 4 


1.0±0.6 


Pi Pen 1 


31 TanD3 


HR 1543 


427.1 


437.8 


0.8 


Clle... 


2904.8 


11. 50 5 


-3.9±0.6 


RTI Anr 


2Q TanD3 


51 Ori 


466.5 


464.3 


0.7 


M8 


3553.5 


3.7 4 


-0.9±0.6 


V A nr 

A Aur 


zyuctuz 


v Aur 


163.8 


166.1 


0.2 


K2 


3664.0 


1.8 4 


2.5±1.1 




23Jan03 


HD 31312 






0.2 




3378.4 




-0.2±0.6 




25Jan03 








0.2 




3429.5 




-0.9±0.6 




28Jan03 








0.2 




3507.5 




-2.1±0.6 


V Mon 


29Jan03 




340.5 


332.7 




M6e 


2771.1 


5.6 4 


2.6±0.6 


W Cnc 


28Jan03 


uj Hya 


393.2 


391.8 


0.3 


M7e 


3415.2 


4.7 4 


-0.7±0.6 


X Hya 


30Jan03 


28 Hya 


301.1 


301.2 


0.6 


M7e 


2451.7 


5.0 4 


-0.2±0.6 


R LMi 


28Jan03 


HR 5464 


372.2 


375.4 


0.3 


M7e 


3530.8 


13.2 2 


-2.6±0.7 


V Boo 


29May02 


p Boo 


258.0 


261.1 


0.6 


M6e 


2100.6 


5.2 4 


-0.7±0.7 




28Jan03 


HR 5464 






0.6 




3448.6 




0.2±0.6 


S CrB 


29May02 


52 Boo 


360.3 


365 


0.7 


M7e 


2121.3 


9.1 2 


0.3±0.8 




28Jan03 


HR 5464 






0.4 




3509.1 


8.9 2 


1.5±0.6 


S Ser 


28Jun02 


k Ser 


371.8 


373.7 


0.9 


M5e 


3557.1 


5.35 6 


-1.3±1.5 


BG Ser 
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k Ser 


143 


386.1 


0.7 


M6me... 


3042.7 


6.71 6 


0.1±0.9 
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CGCVS AAVSO 



Target 


Date 
UT 


Calib. 


Period 

(days) 


Period 


Phase 


Sp. Ty. 

(cm) 


Pmax 

(mas) 


(deg) 








30Jan03 


110 Vir 






0.3 




3328 


.4 




0. 


,7±0. 


.7 
















3305 


,8 




-0 


,2±0. 


.7 


R Ser 


27May02 


k Ser 


356.4 
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0.4 


M7IIIe 


2089. 


.8 


7.6 2 





.0±1 


.0 




28May02 


k Ser 






0.4 




2052 


.9 







,4±0. 


.8 


V CrB 


30Jan03 


HR 5464 


357.6 


361.3 


0.9 


N... 


3271. 


.8 


7.26 5 


-1 


,5±0. 


.6 


RU Her 


25May02 


HR 5947 


484.8 


494.1 


0.5 


M7e... 


2116. 


.5 


8.71 6 


1 


.3±0. 


.9 


RT Oph 


28Jun02 


(3 Oph 


426.3 


424.8 


0.0 


M7 


3500. 


.1 


6.52 6 


-1 


.0±1. 


.0 


X Oph 


26May02 


7 Aql 


328.9 


341.1 


0.4 


K1III+... 
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.1 


12.97 3 


-0 


,6±0. 


.8 


R Aql 


26May02 


7 Aql 


98/1 9 


979 Q 


u.o 


M7TTTc. 

ivi i iiie 
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.2 


Q "X 2 

y.o 


1 
-1. 


1±0. 


.8 


W Aql 


28May02 


7 Aql 


490.4 


480.8 


0.1 


S:... 


1650. 


.0 


11.08 5 


1 


.2±1. 


.0 


RT Aql 


24Jun02 


7 Aql 


327.1 


325.7 


0.0 


M7e 


3488. 


.8 


7.24 6 


1 


.9±0. 


.8 


BG Cyg 


29Jun02 


7 Sge 
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285.4 


0.1 


M7e 


3793 


.5 


4.14 6 


-1 


,8±1. 


.0 


RR Aql 


28May02 


7 Aql 


394.8 


398.7 


1.0 


M7e 


1795. 


.9 


10.73 6 





,2±1. 


.0 


UCyg 


28Jun02 


P Oph 


463.2 


469 


0.0 


R... 


3732 


.1 


7.05 5 


-3 


.5±1. 


.4 


V1426 Cyg 


30Oct02 


P Cyg 


470 


481.2 


0.9 


C 


3816. 


.0 


10.8 7 


-0 


.8±0. 


.9 


R Peg 


31Oct02 


7 Psc 


378.1 


378.2 


0.8 


M7e 


3579. 


.7 


7.0 2 


-1 


.9±1. 


.2 


(b) Semi-regular Variables: 






















p Per 


21Nov02 


24 Per 


50 






M4II 






15.53 8 


1 


.3±3 


.5 


W Ori 


31Oct02 


40 Ori 


212 


(446) 




CII... 


3394 


.3 


9.7 7 


5. 


.2±1. 


.8 


CE Tau 


30Oct02 


HD 33554 
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M2Iab: 


3717. 


.4 


9.1 9 


-4. 
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.2 



Table 2 — Continued 



CGCVS AAVSO 



Target 


Date 
UT 


Calib. 


Period 

(days) 


Period 


Phase 


Sp. Ty. 

(cm) 


(mas) 


9jJ£) 

(deg) 






X Cnc 


19Jan03 


uj Hya 
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379.6 


0.3 


CII... 


3527.5 


7.62 10 


1. 


8±0.6 




25Jan03 








0.3 




3473.1 




-0. 


,4±0.6 




28Jan03 








0.3 




3410.5 




2. 


1±0.6 


T Cnc 


23Jan03 


uj Hya 


428 


499.5 




N... 


3447.1 


6.6 11 


-2. 


0±0.6 




28Jan03 












3409.6 




0. 


3±0.6 


U Hya 


25Jan03 


a Crt 


450 






CII... 


2902.0 


10.8 4 


-2. 


4±0.7 




30Jan03 


HR 5464 










2489.6 




-2. 


.8±0.6 


V Hya 


23Jan03 


a Crt 


530.7 




1.0 


c... 


2640.1 


13.0 4 


2. 


1±0.6 


Y CVn 


29May02 


HR 4690 
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Clab:... 


2022.7 


11.6 7 


-3. 


6±1.0 


RT Vir 


28May02 


v Boo 
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M8III 


1902.3 


12.38 12 


0. 


7±0.7 


SW Vir 


28May02 


v Boo 




1 ^ 4 




M7TTT 

1V1 l 111 




1 fi 94 12 

1U .Zi^i 





s+n q 

. oznu . v 


RX Boo 1 


23May02 


p Boo 


340 






M7.5 


2112.2 


17.48 12 


-3. 


9±1.0 


ST Her 


23May02 


HR 5763 


148 






M6s 


1945.5 


9.3 13 


-0. 


1±1.0 


X Her 


25May02 




95 






M 




12.1 13 






g Her 1 


29May02 


52 Boo 


89.2 






M6III 


2063.0 


12.67 12 


-2. 


5±1.0 


R Lyr 1 


25May02 


HR 6695 


46 






M5III 


2056.4 


13.3 13 


0. 


3±0.9 


V Aql 


28Jun02 


(3 Oph 


353 


377.1 




CII... 


2870.1 


10.1 7 


-0. 


8±1.0 


EU Del 


26May02 


7 Sge 


59.7 




0.7 


M6III 


2063.7 


9.8 9 


-0. 


,8±0.8 


(c) Irregular Variables: 


















Del2 Lyr 
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M4II 


2099.3 


10.32 14 


-0 


.6±0.9 
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CGCVS 


AAVSO 








Target 


Date 


Calib. 


Period 


Period 


Phase Sp. Ty. 


Pmax 






UT 




(days) 




(cm) 


(mas) 


(deg) 


EPS Peg 


25Jun02 


7 Peg 






K2Ib 


3657.3 


7.7 13 -O.lil.l 


TX PSc 1 


25Jun02 


7 Peg 




(255.5) 


CII... 


3571.0 


11. 2 7 4.9±1.6 



■>■>■>■> ■>■>■>■> ■>■>■>■> 



1 At least one of the baseline resolve this target to the level of below 5% in V 2 . 
2 Millan-Gabet et al. (2005) 
3 Thompson et al. (2002) 
4 van Belle (1999) 
5 van Belle et al. (1997) 

6 van Belle, Thompson & Creeck-Eakman (2002) 

7 Dyck, van Belle & Benson (1996) 

8 Di Benedetto (1993) 

9 Dyck, van Belle & Thompson (1998) 

10 Richichi & Calamai (1991) 

11 Richichi, Lisi, Calamai (1991) 



Mennesson et al. (2002) 
Dyck et al. (1996) 
Sudol et al. (2002) 

The uncertainties indicates random errors only (see text). 
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Out of 56 AGB stars, 16 are well resolved (i.e. at least one baseline gives a V 2 mea- 
surement less than 5% of the point source value) from our observations. Among these 16 
targets, 12 of them show asymmetry from our observations. Thus, if we consider only well 
resolved targets, 75% of AGB stars show asymmetry. Targets from our measurements that 
are well resolved are marked in Tables 1 & 2 (see the footnote). Interestingly, all well resolved 
targets that don't show asymmetry (except R Lep), namely, g Her, R Lyr & RX Boo are 
SR variables; R Lep is a carbon Mira. Thus, if we consider only well-resolved oxygen rich 
Mira stars then our asymmetry detection is 100% and in the case of SR variables the suc- 
cess rate is 40%. The well-resolved Mira variable R LMi did not show asymmetry from our 
January 2003 observations. However, we detected asymmetry in this target during March 
2003 observations. Similarly, the well-resolved Irr variable TX Psc did not show asymmetry 
from our June 2002 measurements. However, we detected asymmetry in TX Psc from our 
follow-up observations taken in October 2002. 

At least five targets with detected asymmetry, namely, S CrB, RU Her, R Aql, V Hya 
and Y CVn have earlier measurements taken with relatively shorter baselines and we did 
not detect asymmetry from these measurements. 

In order to understand the role of angular resolution on the asymmetry detection, we 
derived the number of pixel elements (Npj X ) in an imaging sense, defined as the angular 
diameter divided by the angular resolution (X/2B max ) for all targets. Here B max is the 
maximum baseline of our observations listed in Tables 1 & 2 for all our targets. We plotted 
our measured closure phase values against number of pixel elements in Figure (3). This figure 
clearly shows that the positive closure phase cases are those that have pixel elements close 
to or greater than unity, meaning they are well resolved. This suggests that the detected 
asymmetry features are probably on the surface of the stellar disk or visible only in projection 
against the stellar disk (such as might be the case for patchy dust opacity at 1.5 to 2.5 stellar 
radii). 



5. Simple Models 

We have assumed a 2-component brightness distribution model in order to find the 
simplest possible implications of the measured closure phase signal. This model consists of a 
uniform-disk star with an intensity distribution I p (r) and an unresolved secondary component 
(bright spot, companion or dust clump) with an intensity distribution I s (f '— 5r), where 5r 
is the separation vector between the optical centers of the components. The total intensity 
is / = I p (r) + I s (r — Sr). The complex visibility is the Fourier transform of the brightness 
distribution. Thus, the complex visibility of this composite object could be written using 
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the shift theorem for Fourier transforms as 

V(g) = V p (g) + V s (g) e ik ^ 

where k = 2n/\, A is the wavelength of observation and g is the baseline vector B. 
The complex visibility for the baseline Bab could be written as 

V AB = V P AR e l ^B + V s J<t>° ABe ikBAB.8r 



(7) 



(8) 



The visibility phase for the baseline B AB is 



<Pab = arctan 



Im V AB 
Re Vab 



■ arctan 



= arctan 



v ab sH4>ab) + v ab sil H<? AB 
Vab cos(0^ b ) + VXb cos(05, 



Vab sin(fe) + V S AB ^Wab + kB AB Sr) 
V'ab cos((f) p AB ) + VX B cos((f) s AB + kB AB -Sr) 

cos(kB A B-^r) + V; 



(9) 



(10) 



) smikBAB-^f) 



(11) 



, COS(0 

/t/i. 

i> s AB ) cos(kB A B-Sr) - \/| B sin(0^ B ) sin(/c5 j4B .5r) 

The individual components of the brightness distributions are assumed to be circularly 
symmetric. Hence, sm((p AB ) = sin(0^ B ) = 0, but cos((fi AB ) and cos((f) AB ) can be +1 or 
-1, depending on details of the case. Thus, 



VXb cos((f) AB )sm(kB AB .5r) 



4>ab = arctan 

'Vab cos ( < I ) ab) + Vab cos((f) AB )cos(kB AB .8r) 
Unresolved Secondary component: 

For an unresolved (i.e. a point source) (f> AB = 0. Thus, 

VXb sm(kB AB .5r) 



(12) 



4>ab = arctan 



V V AB cos(0%) + VXb cos(kB AB .Sr) 



(13) 



For a uniform-disk star, 



VL 



2J 1 (7riVff/2) 



1 P) 



(14) 
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and 



Kb cos((p p AB ) 



(15) 



For an unresolved secondary component, 

Vab = Is, 



(16) 



where I p and I s are the normalized star and secondary component intensities (i.e. I p + I s 
= 1), and 



N AB = &UD 

ptx (X/2B AB ) 



(17) 



Combining eqn (13), (15) and (16), we get 



4>ab = arctan 



I s sm(kB AB .5r) 



I P ( 2Jl ^|f 2 /2) ) + /. cos(kB AB Sr) 



(18) 



Resolved Secondary component: 

For a resolved Gaussian secondary component (such as dust clump) (p AB = 0, and 

Vis = [exp(-(n(3N^/2a) 2 )] I s (19) 
where, a = 2 yjln{2) and f3 is the ratio of the size of the primary component to the secondary. 
Thus, 



*>ab = arctan 



I a exp(-(7r(3N^/2af) sm(kB AB Sr) 
l„ ( 2Jli : N A T, m ] + I s exp{-{-K(3N^/2af) cos(kB AB Sr) 



(20) 



For a resolved uniform-disk secondary component (such as stellar companion) 4> S AB = 
or ±7r, and 



V s -- 

V AB — 



23^^/2) 



Is, 



and 



V S AB cos(^) 



2J x WNgl2j 



Is, 



(21) 



(22) 
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Now, 



4>ab = arctan 




sin(kBAB-Sr) 



— * — * 

cos(kB AB .8r) 



(23) 



Using eqn (4), (18), (20) and (23), we generate closure phase models for the IOTA 
configuration A35B15C00 (i.e. Be a = 35 hi; Bbc — 15 m) which are shown with our 
data in Figure 3. We explored the parametric space with four cases, namely (1) surface- 
unresolved spot, (2) surface-unresolved companion, (3) surface-resolved companion and (4) 
surface-resolved Gaussian dust clump. The values assumed for various parameters of these 
models are primarily for the purpose of illustration. More thorough treatment of the physical 
parameters chosen in the models will be presented in the forthcoming article. We have shown 
two models for each of the four cases (totally eight models) in Figure 3. The two models 
differ only in the position angle of the secondary feature - one assumes the secondary feature 
at the position angle of 21.8° and the other at the position angle of 201.8° (i.e. a 180° 
rotation). The reason for choosing this axis (along the direction of the largest baseline 
of the IOTA array) is that the closure-phase signal is approximately maximum when the 
secondary feature is assumed along this axis. The brightness distribution of the primary 
component is assumed to be an uniformly illuminated disk. The flux of the secondary 
component is assumed to be 3% of the total flux for surface-unresolved secondary cases 
and 30% for surface-resolved secondary component cases; the reason for choosing these flux 
values for the companion is that the corresponding models compare well with our closure 
phase measurements. The spot models (case 1) assume a surface-unresolved bright spot 
at the edge of the stellar disc (5r = 8jjd/2)- The surface-unresolved binary models (case 
2) assume a surface-unresolved companion or dust clump at 5 stellar radii (8r = 2.5 6 UD ). 
The surface-resolved binary models (case 3) assume companion with (3 = 0.99 (UD angular 
size of the secondary component is 99% of the primary) at 5 stellar radii (5r = 2.5 6ud)- If 
companion stars provide enough light to produce detectable asymmetry near maximum light, 
then they should also produce a wider, flattened, minima in the light curves (Merrill 1956). 
While at least two of the stars are known binaries and relatively mild symbiotic systems (R 
Aqr and o Ceti) none of the stars show a filled-in minimum on the AAVSO light curves. The 
dust models (case 4) assume a Gaussian shaped dust clump ((3 = 1; same equivalent size as 
the primary component) at 5 stellar radii (5r = 2.5 Quo)- As can be seen in Figure 3, the 
unresolved spot models compare well with the observed data and the unresolved companion 
models are not as good as the unresolved spot models in explaining the observed closure 
phase data. The resolved secondary feature models (case 3 & 4) show that the secondary 
components have to be significantly brighter (may be physically unreasonable) in order to 
produce detectable closure phase signals and even so, the fits at low N pix are poor. 
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Fig. 3. — Measured closure phases are plotted against number of pixel elements (see the text). Targets 
with positive asymmetry detection are shown as diamond symbols and non-detection measurements are 
shown as plus symbols. Targets that are well resolved frequently show large closure phase. The solid lines 
refer models where the secondary feature is assumed at the position angle of 21.8° (along the direction of 
the largest baseline of the IOTA array) and the dotted lines refer to models where the secondary feature is 
assumed at the position angle of 201.8° (21.8° + 180°). The flux of the unresolved secondary components are 
assumed to be 3% and that of the resolved secondary component is assumed to be 30%. Top left: Unresolved 
spot at the edge of the stellar disc (Sr = Oud/^)] Top right: Unresolved companion at 5 stellar radii (Sr 
= 2.5 6ud)', Bottom left: Resolved UD companion at 5 stellar radii (Sr — 2.5 Bud)] The diameter of the 
secondary is assumed to be 99% of the diameter of the primary; Bottom right: Resolved Gaussian dust 
clump at 5 stellar radii (Sr — 2.5 9ud)', The equivalent size of the secondary is assumed to be same as that 
of the primary. 
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6. Discussion 

The significance of the present results (i.e., that 1/3 are asymmetric and 2/3 are not) 
depends on what causes the asymmetry, and that is not yet known. As noted in the intro- 
duction, although planetary nebulae are predominantly bipolar, a large fraction (possibly 
even all) of proto-PN are axisymmetric inside an apparently spherical AGB wind remnant, 
suggesting that the PN asymmetry has arisen only after or as the star left the AGB. Asym- 
metries have been reported before for Miras, but these are mostly for isolated examples, or 
for maser emission that is very sensitive to the local conditions and thus will tend to exag- 
gerate any physical departure from spherically symmetric flow. The present results, referring 
to intensity near the stellar flux maximum, are much less sensitive to small variations in the 
conditions. It's worth mentioning here that the non-detections don't preclude asymmetries. 
It could be that they are just not resolved or the asymmetries are too small. 

The SiO maser emission also arises well above the photosphere (e.g. Humphreys et al 
2002; Cotton et al. (2004)); in fact, there is, to our knowledge, no report yet of asymmetry 
that can be assigned unambiguously to the stellar surface. With multiple narrow-band mea- 
surements at carefully selected wavelengths, Perrin et al (2004) have shown that it is possible 
to disentangle photospheric and shell contributions. The results reported here suggest that 
we will also be able to sort out some shape information in the next generation of observations 
as well as separate the photospheric contribution from the circumstellar one. 

There is a considerable literature concerning the non-circular and non-spherical sym- 
metry common among planetary nebulae; see, for example, the proceedings of Asymmetric 
Planetary Nebulae III (2004, ASP Conf. Ser. 313). Mechanisms may be roughly divided 
into deep and superficial. Deep mechanisms include internal convective structure with large 
convection cells (proposed by Schwarzschild (1975) on the basis of simple scaling arguments), 
non-radial pulsation, and/or rotation. 

These stars have massive envelopes and large radii; no reasonable reservoir of angular 
momentum other than incorporation of a relatively massive (> 0.1 M ) companion will pro- 
vide sufficient angular momentum for rotational asymmetry or the usual non-radial-pulsation 
associated with rotation. Large convective cells might stimulate non-radial modes in the ab- 
sence of significant rotation, or might lead to modulations in the surface brightness from 
rising or falling elements. Evolutionary models show a radiative layer above the convective 
layer in these stars (Ostlie & Cox 1986) and the scale heights at the photosphere are only 
a percent or so of R*; these facts suggest that convective modulation will be on a smaller 
scale and with less contrast than is needed to explain these observations, but more detailed 
modeling should be done before the possibility of convection-based modulation is ruled out. 
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Superficial or atmospheric mechanisms include magnetic structures (e.g. Soker & Zoabi 
2002; Blackman et al. 2001), discrete dust cloud formation as for R CrB stars, and inter- 
action of a planet or companion with the stellar wind (Struck, Cohanim & Willson 2004; 
Mastrodemos & Morris 1998, 1999). The conclusion that perhaps all of the Miras show some 
asymmetry while only about half of the non-Miras do may be understood in a couple of ways. 
Miras comprise a well-defined subset of long period variables, namely those with large visual 
amplitudes, relatively regular variation, cool effective temperatures, and moderate progeni- 
tor masses. SR classes differ in visual amplitude (SRa), degree of irregular variation (SRb), 
warmer effective temperature (SRd), and higher progenitor mass (SRc). Within each of 
these classes there are further probable subclasses. Our current understanding is that the 
high visual amplitude is partly the result of variable atmospheric opacity (Reid & Goldston 
2002); this variable opacity is closely tied to the fact that these stars are losing mass at a 
high rate (>10~ 7 M Q per year up to 10 5 M per year) so the evolutionary status of Miras 
is that they are stars entering the final "superwind" (massive outflow) stage on the AGB 
(e.g. review Willson 2000). A number of stars initially classified as Miras are reclassified as 
SRb when their light curves develop irregularities, making the boundary between these two 
classes somewhat fuzzy. Similarly, most relatively regular carbon star LPVs have smaller 
amplitudes in the visual than do the oxygen-rich stars, and this may be telling us more about 
the sources of atmospheric opacity than about the evolutionary state, so that the SRa-Mira 
boundary is also fuzzy. Thus one interpretation would be that the asymmetry shows up 
when there is a sufficiently massive outflow to produce the large Mira amplitude for the 
oxygen-rich stars, and that the SR variables with asymmetries are those with different visual 
opacity but similarly massive outflows. Either the same mechanism leads to outflow and 
asymmetry (e.g. non-radial pulsation) or the outflow sets up conditions for asymmetry to be 
seen. In the first case, the non-radial structure originates at the photosphere; in the second, 
with aperiodicities in the outflow. 

Most Miras are surrounded by translucent "molecular shells", a locus in the outflow 
where molecules and probably dust provide high local opacity, whose IR and visual optical 
depth is on the order of 1 (0.1 to several) - Perrin et al (2004). The physics of dust formation 
in the context of large-amplitude pulsation and consequent shocks is reviewed in Willson 
(2000). Dust grains nucleated in the refrigerated zone between shocks may require several 
pulsation cycles to grow to sufficient size to generate an outflow, and this would naturally 
lead to critical dust levels appearing in different cycles at different positions around the 
star. Whether by this or another mechanism, the translucent shell is likely to have a patchy 
opacity, allowing more of the photospheric light through in some places than in others. Thus 
a plausible explanation for the possibly universal asymmetry in Miras would be the formation 
of an inhomogeneous translucent molecular screen around 1.5 to 2.5 stellar radii. 
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In conclusion, we carried out a survey of AGB stars with the IOTA 3-telescope imaging 
interferometer at near-infrared wavelengths, searching for asymmetry in their flux distribu- 
tions. We find that 29% of our sample show asymmetry. If we restrict the sample to only 
well resolved targets, then 75% of AGB stars, 100% of oxygen-rich Mira stars show asym- 
metry from our observations. On this basis, we hypothesize that all Mira stars might show 
detectable asymmetry if observed with adequate spatial resolution. The large frequency of 
asymmetry reported here suggests that angular size measurements and limb darkening stud- 
ies of AGB stars carried out with 2-telescope optical long baseline interferometers should 
be interpreted with caution. We have initiated a systematic mapping program, namely, 
'The Mira Imaging Project' funded by National Science Foundation (NSF) at the IOTA, ISI 
and VLBA interferometers to connect the asymmetry in space and time, and pin-point the 
mechanism(s) responsible for observed asymmetry. 
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